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Abstract: The Time Domain Integral Equation (TDIE) method is an appealing method for time- 
dependent scattering analysis of electromagnetic waves. As a surface integral equation method, it 
is computationally efficient for scattering models at arbitrarily shaped objects and automatically 
satisfies radiation conditions for exterior problems. The time-domain formulation allows for the 
simulation of broadband signals at materials with nonlinear constitutive relations. However, its 
applicability has long been restricted by computational difficulties to obtain stable simulations for 
large-scale objects. Here, a stability analysis will be presented that is based on an established 
functional framework and stability theorem for a variational formulation of the electric field 
integral equation (EFIE). The functional framework will be extended to the differentiated version 
of the EFIE, which is more efficient to discretize and more conventional in engineering literature. 
Then, a discrete equivalence between Petrov-Galerkin and collocation schemes will be presented. 
This allows for a stability study of popular temporal basis functions. The impact of the stability 
theorem on the choice of discretization scheme is given by the recommendation to use quadratic 
spline basis function in the standard Marching-on-in-Time scheme. Computational experiments of 
time-dependent scattering at a generic aircraft confirm the stability analysis of the TDIE method. 
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TDIE simulation, excitation 





| «a ZA yy: 


NE 


v n ~ ds i? 


n i U a W 
i HU | p LE n i WW "M A T TRE 
ra L ze ij mj 



















NLR - Dedicated to innovation in aerospace 2 


TDIE simulation, electric surface current and 
scattered field 
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Formulation 


u OJ(r,T) ão V -J(r, t)dt ia 
manx || (R = — UV z dr = E'(r,t) (1) 


u O*J(r ,V'J(r,T) 
-nxmx q || (EE E 2) dr = 


J(r,t): unknown electric surface current 
I’: scattering surface 
Tr = t — R/c: retarded time 

= |R| = |r — r'| 


(1): original EFIE 
(2): differentiated EFIE 


OE"(r’,t) 


at (2) 
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Boundary element method 





Discrete variational formulation of £(J) — EM: 


(EE: Jig) dh) = (E" hj) Vj € {1,2,..., N}. 
Space-time approach: 


g,(r,t) = f„(rju,(t), 
h;(r, t) = fa(r)ve(t), 


with f(r) RWG functions. DES 
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Marching-on-in-Time scheme 


Discrete time axis: t; = 9 At. 


Test functions: v(t) = (t — At) (collocation). 
Basis functions: u(t) = T(t — kAt) with 


Fo(s), —At < s < 0, 
Fls), O<s< At, 


Fals), (d—1)At < s < dAi, 
O, else, 


plecewise polynomials of degree d. 
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Marching-on-in-Time scheme 





Zol; = V; — SZ, fora=1,2,...,N; 
j=1 
MoT scheme for differentiated EFIE: 
e causality, 
e accuracy, 
e efficiency. 
Challenge: stability. 
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Stability theorem, current status 





The variational formulation of the original EFIE 
has a unique, bounded solution for 


Ftest — H-?, 
F basis — H? 
with 
H = HS (Ra; H-3 (div, r) | 


[1] |. Terrasse, “Résolution mathématique et numérique des équations de Maxwell instationnaires 
par une méthode de potentiels retardés,” Ph.D. dissertation, École Polytechnique, 1993. 
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Stability theorem, example functions 





LALA 


Dirac delta: HT >. Step: H 2. Hat: H2. Quadratic spline: H 2. 
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Stability theorem, extensions 





Required extensions for use in MoT schemes: 
1. Continuous level: original/differentiated EFIE 
2. Discrete level: Dirac delta as test function 
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Stability theorem, model extension 


original EFIE — differentiated EFIE 


Variational formulation of original and differentiated EFIE: 


(8, Eorig(J)), = (g, ES. Vg E Forig,test, 
(h, Eaitt(J)) = (h, SES Vhe F diff test. 


— Oh A 
We have proven that for g = + — 20h: 
e the variational formulations are equivalent and 


e g EH eheh. 
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Stability theorem, scheme extension 


space-time Galerkin — MoT 


Discrete model equation: 





Zol; = Vi — Zili-1 ala — +++ — Zpli-p 


with discretization matrices Z; and discrete solution I; at 
time t;. 


If ay” = Zo then = m | kad 
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Stability theorem, scheme extension 


Examples of equivalencies: 


MoT test MoT basis stG test stG basis 


H, 
H, 
ux 


D: 


+ 


SIE 
ple 


A 


hat spline Lagrange 


o 
E 
E 


Dirac 
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Stability theorem, summary 


Stability of the differentiated EFIE: 


Ftest — H~? and J basis = H2. 
Example: 


fstep E Ftest aNd fhat E Fbasis. 
Discretely equivalent with: 


fDirac £ J test and Íspline E F| basis - 


EN] EN 


Dirac step hat spline 
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Nuance: quadrature accuracy 








Recommendations 


Use: 

. the differentiated EFIE, 

. RWG functions, 

. the MoT scheme, 

. Quadratic spline basis functions, 

. Quasi-exact integration, 

. an oversampling as time step size. 


2 
3 
4 
D 
6 


NLR - Dedicated to innovation in aerospace 





Bistatic RCS of an aircraft 








MoM 
-------- TDIE 
Aircraft of size 60 m. 
Triangular surface mesh 
with 6756 edges. A azimuth 


À = 10m 
+x-propagation 
+z-polarization 


180 
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Conclusion 


The MoT scheme with quadratic spline basis 
functions fit within the stability theorem of 
space-time Galerkin schemes. 
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